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ABSTRACT 

We measure the angular power spectrum Cg of radio galaxies in the NRAO VLA Sky 
Survey (NVSS) using two independent methods: direct spherical harmonic analysis 
and maximum likelihood estimation. The results are consistent and can be understood 
using models for the spatial matter power spectrum and for the redshift distribution 
of radio galaxies at mJy flux-density levels. A good fit to the angular power spec- 
trum can only be achieved if radio galaxies possess high bias with respect to mass 
fluctuations; by marginalizing over the other parameters of the model we derive a 
68% confidence interval 1.53 < 6o x er 8 < 1.87, where bo is the linear bias factor for 
radio galaxies and describes the normalization of the matter power spectrum. Our 
models indicate that the majority of the signal in the NVSS Ci spectrum is generated 
at low redshifts z<0.1. Individual redshifts for the NVSS sources are thus required 
to alleviate projection effects and probe directly the matter power spectrum on large 
scales. 
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, I 1 INTRODUCTION 
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Active Galactic Nuclei (AGN) mapped in radio waves are an 
interesting probe of large-scale structure. They can be rou- 
tinely detected out to very large redshift (z ~ 4) over wide 
areas of the sky and hence delineate the largest structures 
and their evolution over cosmic epoch. Radio emission is in- 
sensitive to dust obscuration and radio AGN are effective 
tracers of mass: they are uniformly hosted by massive ellip- 
tical galaxies and have been shown to trace both clusters 
(Hill & Lilly 1991) and superclusters (Brand et al. 2003). 

The current generation of wide-area radio surveys such 
as Faint Images of the Radio Sky at Twenty centimetres 
(FIRST; Becker, White & Helfand 1995) and the NRAO 
VLA Sky Survey (NVSS; Condon et al. 1998) contain radio 
galaxies in very large numbers (~ 10 6 ) and have allowed 
accurate measurements of the imprint of radio galaxy an- 
gular clustering. These patterns are considerably harder to 
detect in radio waves than in optical light due to the huge 
redshift range that is probed. Whilst this provides access 
to clustering on the largest scales, it also washes out much 
of the angular clustering signal through the superposition 
of unrelated redshift slices. The angular correlation function 
was measured for FIRST by Cress et al. (1996) and Maglioc- 
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chetti et al. (1998) and for NVSS by Blake & Wall (2002a) 
and Overzier et al. (2003). 

The NVSS radio survey, covering ~ 80 per cent of the 
sky, permits the measurement of fluctuations over very large 
angles. Blake & Wall (2002b) detected the imprint of the 
cosmological velocity dipole in the NVSS surface density, 
in a direction consistent with the Cosmic Microwave Back- 
ground (CMB) dipole. In this study we measure the angular 
power spectrum, Ci, of the radio galaxy distribution (Balei- 
sis et al. 1998). This statistic represents the source surface- 
density field as a sum of sinusoidal angular density fluctua- 
tions of different wavelengths, using the spherical harmonic 
functions. The angular power spectrum is sensitive to large- 
angle fluctuations and hence complements the measurement 
of the angular correlation function, w(8), at small angles. 

Measurement of the Ce spectrum has some advantages 
in comparison with w(8). Firstly, the error matrix describing 
correlations between multipoles £ has a very simple struc- 
ture, becoming diagonal for a complete sky. This is not the 
case for the separation bins in a measurement of w(9): even 
for a full sky, an individual galaxy appears in many separa- 
tion bins, automatically inducing correlations between those 
bins. Secondly, there is a natural relation between the an- 
gular power spectrum and the spatial power spectrum of 
density fluctuations, P(k). This latter quantity provides a 
very convenient means of describing structure in the Uni- 
verse for a number of reasons. Firstly its primordial form 
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Figure 1. NVSS sources with S1.4GHZ > 200 mjy in an equal-area projection. The Galactic plane and Galactic latitudes ±5° are also 
plotted; sources within this region are masked from our large-scale structure analysis as many are Galactic in origin. 



is produced by models of inflation, which prescribe the ini- 
tial pattern of density fluctuations 5p/p. Furthermore, in 
linear theory for the growth of perturbations, fluctuations 
described by different wavenumbers k evolve independently, 
enabling the model power spectrum to be easily scaled with 
redshift. The physics of linear perturbations are hence more 
naturally described in Fourier space. 

In contrast, the angular correlation function is more 
easily related to the spatial correlation function £(r), the 
Fourier transform of P(k). Correlation functions more nat- 
urally serve to describe the real-space profile of collapsing 
structures evolving out of the linear regime. We emphasize 
that although the two functions Ce and w(8) are theoret- 
ically equivalent - linked by a Legendre transform - this 
is not true in an observational sense. For example, w(8) can 
only be successfully measured for angles up to a few degrees, 
but Ci depends on w(9) at all angles. 

We derive the angular power spectrum using two inde- 
pendent methods. Firstly we apply a direct spherical har- 
monic estimator following Peebles (1973). Secondly, we use 
maximum likelihood estimation, commonly employed for de- 
riving the angular power spectra of the CMB temperature 
and polarization maps. These two methods are described 
in Section 3. We find that these two approaches yield very 
similar results (Section 4), which is unsurprising given the 
wide sky coverage of the NVSS. In Section 5 we interpret the 
NVSS angular power spectrum in terms of the underlying 
spatial power spectrum of mass fluctuations and the radial 
distribution of radio sources. Finally in Section 6 we employ 
these models to derive the linear bias factor of NVSS radio 
galaxies by marginalizing over the other model parameters. 



2 THE NVSS RADIO SURVEY 

The 1.4 GHz NRAO VLA Sky Survey (NVSS; Condon et al. 
1998) was performed at the Very Large Array over the period 
1993 to 1996 and covers the sky north of declination —40°. 
The source catalogue contains r; 1.8 x 10 6 entries and is 99 
per cent complete at integrated flux density Si. 4 ghz = 3.5 
mjy. The full width at half-maximum of the synthesized 



beam is 45 arcsec; the majority of radio sources are thus 
unresolved. The relatively broad NVSS beam yields excellent 
surface brightness sensitivity and photometric completeness. 

Before analyzing the survey for large-scale structure we 
imposed various angular masks. Firstly we excluded cata- 
logue entries within 5° of the Galactic plane, many of which 
are Galactic in origin (mostly supernova remnants and HII 
regions). The contribution of foreground Galactic sources at 
latitudes |6| > 5° is negligible. We also placed 22 masks 
around bright local extended radio galaxies contributing 
large numbers of catalogue entries, as described in Blake 
& Wall (2002a). Figure 1 plots the remaining sources with 
5*1.4 GHz > 200 mjy. The remaining NVSS geometry corre- 
sponds to 75 per cent of the celestial sphere. 

Blake & Wall (2002a) demonstrated that the NVSS suf- 
fers from systematic gradients in surface density at flux- 
density levels at which it is complete. These gradients - 
corresponding to a ~ 5 per cent variation in surface den- 
sity at a threshold of 3 mjy - are entirely unimportant for 
the vast majority of applications of this catalogue. How- 
ever, they have a significant influence on the faint imprint 
of large-scale structure. If left uncorrected a distortion of the 
measured angular power spectrum would result, because the 
harmonic coefficients would need to reproduce the system- 
atic gradients as well as the fluctuations due to clustering. 

Blake & Wall (2002a) found that these surface gradients 
are only significant at fluxes Siagh* < 10 mjy. At brighter 
fluxes the survey is uniform to better than 1 per cent - suf- 
ficient to allow the detection of the anticipated cosmological 
velocity dipole (Blake & Wall 2002b). Hence we simply re- 
stricted our Ci analysis to fluxes Si.4GHz > 10 mjy. The 
NVSS source surface density at this threshold is ctq = 16.9 
deg" 2 . 

We note that the broadness of the radio luminosity func- 
tion ensures that the projected clustering properties of radio 
galaxies are not a strong function of flux density in the range 
3 mjy < Si. 4 ghz < 50 mjy, as verified by the correlation 
function analyses of Blake & Wall (2002a) and Overzier et 
al. (2003). In this flux-density range, the redshift distribu- 
tion of the radio galaxies does not change significantly. 
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3 ESTIMATING THE ANGULAR POWER 
SPECTRUM: METHODS 

3.1 Definition of the angular power spectrum 

A distribution of galaxies on the sky can be generated in two 
statistical steps. Firstly, a density field a(8, <f>) is created; 
this may be described in terms of its spherical harmonic 
coefficients ae m : 



(1) 



where Ye m are the usual spherical harmonic functions. Sec- 
ondly, galaxy positions are generated in a Poisson process 
as a (possibly biased) realization of this density field. 

The angular power spectrum Ce prescribes the spherical 
harmonic coefficients in the first step of this model. It is 
defined over many realizations of the density field by 

<\a tm \ 2 >=C e (2) 

The assumption of isotropy ensures that < \a lm \ 2 > is a 
function of only £, not m. 

The angular power spectrum Ce is theoretically equiva- 
lent to the angular correlation function w(8) as a description 
of the galaxy distribution. The two quantities are connected 
by the well-known relation 



Ce = 2ira 



w(6)P e (cose) d(cos6) 



(3) 



where ao is the source surface density and Pe is the Legendre 
polynomial. However, the angular scales on which the signal- 
to- noise is highest are very different for each statistic. w(0) 
can only be measured accurately at small angles up to a few 
degrees (Blake & Wall 2002a), beyond which Poisson noise 
dominates. By contrast, Ce for galaxies has highest signal- 
to-noise at small t, corresponding to large angular scales 
6 ~ 180° /£. Hence the two statistics are complementary, 
the Ce spectrum probing fluctuations on the largest angular 
scales. 



3.2 Spherical harmonic estimation of Ce 

Peebles (1973) presented the formalism of spherical har- 
monic analysis of a galaxy distribution over an incomplete 
sky (for refinements see e.g. Wright et al. 1994; Wandelt, 
Hivon & Gorski 2000). For a partial sky, a spherical har- 
monic analysis is hindered by the fact that the spherical 
harmonics are not an orthonormal basis, which causes the 
measured coefficients ae m to be statistically correlated, en- 
tangling different multipoles of the underlying Ce spectrum. 
However, for the case of a survey covering ~ 80% of the 
sky, the repercussions (discussed below) are fairly negligible, 
implying shifts and correlations in the derived power spec- 
trum that are far smaller than the error bars. We employed 
the original method of Peebles (1973) with only one small 
correction for sample variance. In Section 3.4 we compare 
spherical harmonic analysis with the technique of maximum 
likelihood estimation and find that the two methods yield 
results in good agreement. 

The spherical harmonic coefficients of the density field 
may be estimated by summing over the N galaxy positions 



Aer, 
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For an incomplete sky, these values need to be corrected for 
the unsurveyed regions, so that an estimate of Ce is 

2 

- - ao (5) 
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(Peebles 1973 equation 50) where ao = N/ Af2 and 

-/ 

-j 



Je 



\Ye m \UQ 



(6) 



(7) 



where the integrals are over the survey area Afi, and are 
determined in our analysis by numerical integration. The 
final term in equation 5 corrects for the finite number of 
discrete sources: for a full sky (Ie m = 0, Jem = 1) we expect 
< \At m \ 2 >= ao in the absence of clustering, i.e. this is the 
power spectrum of the shot noise. Note that ao is the appar- 
ent source density N/ AO, not the average over an imagined 
ensemble of catalogues. 

We determined the angular power spectrum for the £th 
multipole, Cj? bs , by averaging equation 5 over m. Because 
the density field is real rather than complex, Cf^ = Ce h l m , 
resulting in I + 1 independent measurements of Cf. 



robs 



e + i 



(8) 



There is no need for us to use the modified weighting formula 
of Peebles (1973 equation 53). In our case, Jem does not vary 
significantly with m. We verified that the modified weighting 
formula produced indistinguishable results. 

One consequence of the partial sky is to "mix" the har- 
monic coefficients such that the measured angular power 
spectrum at £ depends on a range of CV around £' = i: 



< C° hs >= y ' Ce' Re 



(9) 



The angled brackets refer to an imagined averaging over 
many realizations of density fields generated by Ce, in ac- 
cordance with equation 2. Peebles showed that y\, Ree' = 1; 
i.e. mixing does not spuriously enhance the measured power 
(this is accomplished by the factor J lm in equation 5). 

For a complete sky, Ree' = See 1 , where S mn = 1 (m = n) 
or (m 7^ n). For a partial sky, the matrix Rev can be 
computed from the geometry of the surveyed region (Hauser 
& Peebles 1973). Figure 2 illustrates the result for the NVSS 
for t — 10 (computed using Hauser & Peebles 1973 equation 
12). The NVSS covers a sufficiently large fraction of the 
sky (75 per cent) that mixing only occurs at the ~ 15 per 
cent level and can be neglected because the underlying Ce 
spectrum is smooth: 



< Cf s >c 



Ce 



We checked the NVSS Ree' matrix for other multipoles I and 
found very similar results. This argument ensures that the 
measured multipoles are statistically independent to a good 
approximation. 

The statistical error on the estimator of equation 5 is 



(10) 
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Figure 2. The range of I probed by a measurement of C£ bs at I = 
10. The incomplete sky means that C°q s is not just sensitive to 
Cio, but also to the surrounding CV's. However, the large fraction 
of sky covered means that this mixing is small. 



(Peebles 1973 equation 81). There are two components of 
the error: 

• "Shot noise" (cto) because the number of discrete ob- 
jects is finite and therefore does not perfectly describe the 
underlying density field. 

• "Cosmic variance" (Ce) because even with perfect sam- 
pling of the density field, there are only a finite number of 
harmonics associated with the £th multipole. 

The error for the m = case in equation 10 is greater 
because Aeo is purely real, rather than complex. In the lat- 
ter case, we are averaging over the real and imaginary parts 
of Aim, two independent estimates of Ce, which reduces the 
overall statistical error by a factor \/2. For a partial sky 
equation 10 is an approximation, because the variance of 
multipoles of given £ depends on the underlying power spec- 
trum at £' 7^ £. As discussed above, this effect is negligible 
for the NVSS. 

The averaging over m (equation 8) decreases the er- 
ror in the observation. Combining the errors of equation 10, 
assuming estimates at different m are statistically indepen- 
dent, the resulting error in C g is 



a(C2 



(cro + C e ) 



VT+2 
£+1 



(11) 



We used Monte Carlo simulations to verify that equation 11 
produced results within 5 per cent of the true error for all 
relevant multipoles. 

Our only addition to the formalism of Peebles (1973) 
was to increase the total variance on the estimate of Ci by 
a factor l// s k y , where / s k y = AQ/4-k is the fraction of sky 
covered (i.e. multiply equation 11 by /sky). This cor- 
rection factor was motivated by Scott, Srednicki and White 
(1994) as a fundamental property of sample variance for a 
partial sky, and is part of the standard CMB formalism (e.g. 
Bond, Efstathiou & Tegmark 1997). For the NVSS geome- 
try, / s ky = 0.75, thus this correction corresponds to a ~ 10 
per cent increase in the error. 



3.3 Correction for multiple-component sources 

Radio sources have complex morphologies and large linear 
sizes (up to and exceeding 1 Mpc) . A radio-source catalogue 
such as the NVSS will contain entries which are different 
components of the same galaxy (for example, the two radio 
lobes of a "classical double" radio galaxy). The broad angu- 
lar resolution of the NVSS beam leaves over 90 per cent of 
radio sources unresolved; however, the remaining multiple- 
component sources have a small but measurable effect on 
the angular power spectrum. 

It is relatively simple to model the effect of multiple- 
component sources on the estimator for Ce described in Sec- 
tion 3.2. The relevant angular scales (£ < 100) are much 
bigger than any component separation, and equation 4 can 
be replaced by 



A, 



i=l 



CiYlm 



where iV ga i is the total number of galaxies and a is the 
number of components of the ith galaxy. Thus the quantity 

< A tm >= Agai c < Y tm >= N < Y lm > 

is unchanged by the presence of multiple components (c de- 
notes the average number of components per galaxy, and 
TV = c x /Vg a i is the total number of catalogue entries, as in 
equation 4). But C£ OT S in equation 5 depends on 

< \ A lm \ 2 > = 7V ga ic2 < \Y lm \ 2 > +N 2 al c 2 < Y im >< Yi m > 

— — N < \Y lm \ > +N < Y tm >< Yl m > 
c 

Multiple-component sources only affect the first term in this 
expression, producing an offset in the Ce spectrum 



ACe = 



A(< \A lm \ 2 >) 



- 1 



N < \Y lm \ 2 > 
Jem 



But Je m =< IVeml > from equation 7, and this expres- 
sion simplifies to an offset independent of £: 

ACe = (i-l Wo 



Most multiple-component sources in the NVSS catalogue 
are double radio sources. Let a fraction e < 1 of the radio 
galaxies be doubles. Then c = 1 + e and c 2 = 1 + 3e, thus 
the constant offset may be written 



ACe « 2 e a 



(12) 



We can deduce e = 0.07 ± 0.005 from the form of the NVSS 
angular correlation function w(8) at small angles 6 < 0.1°, 
where double sources dominate the close pairs (see Blake 
& Wall 2002a and also Section 4.2). This correction was 
applied to the measured NVSS Ce spectrum and successfully 
removed the small systematic offset in Ce at high I. 

3.4 Maximum likelihood estimation of Ce 

A sophisticated suite of analytical tools has been developed 
by the CMB community for deriving the angular power 
spectra of the observed CMB temperature and polariza- 
tion maps. These methods can also be exploited to ana- 
lyze galaxy data (see for example Efstathiou & Moody 2001, 
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Figure 3. Maximum likelihood estimation assumes a Gaussian 
distribution for the pixclizcd data. This histogram shows the dis- 
tribution of values of Xi (see equation 13) for the NVSS 10 mjy 
sample, pixelized using HEALPIX parameter n s id e = 32. The 
overplotted Gaussian curve demonstrates that Gaussianity is a 
good approximation for this sample. 



Huterer, Knox & Nichol 2001 and Tegmark et al. 2002). In 
this approach the power spectrum is determined using an 
iterative maximum likelihood analysis, in contrast to the di- 
rect estimator discussed in Section 3.2. The likelihood is a 
fundamental statistical quantity, and this analysis method 
permits straightforward control of such issues as edge effects, 
noise correlations and systematic errors. 

The starting point for maximum likelihood estimation 
(MLE) is Bayes' theorem 

P{a\DI) oc P{a\I) P{D\aI) 

where a are the parameters one is trying determine, D is 
the data and I is the additional information describing the 
problem. The quantity P{a\I) is the likelihood, i.e. the prob- 
ability of the data given a specific set of parameters, while 
the left-hand side is the posterior, i.e. the probability of the 
parameters given the data. 

We will assume that the sky is a realization of a sta- 
tionary Gaussian process, with an angular power spec- 
trum Ce- We assume no cosmological information about 
the distribution of the Cg. The rendition of the sky will 
be a pixelized map, created by binning the galaxy data 
in equal-area cells such that the count in the ith cell 
is m, effectively constructing a "temperature map" of 
galaxy surface density. We performed this task using the 
HEALPIX software package (Gorksi, Hivon & Wandelt 
1999; http://www.eso.org/science/healpix). We chose 
the HEALPIX pixelization scheme n s id c = 32, which cor- 
responds to 12,288 pixels over a full sky. The angular power 
spectrum may be safely extracted to multipole ^ ma x ~ 
2 x n s id c . We then defined a data vector: 



- V± — i 

n 



(13) 



where n is the mean count per pixel. Figure 3 demon- 
strates that the data vector Xi for the NVSS sample is well- 
approximated by a Gaussian distribution, as assumed in a 
maximum likelihood analysis. 

The covariance matrix Cfj due to primordial fluctua- 
tions is given by 



where Pe is the Legendre polynomial and Oij is the angle be- 
tween pixel pair {i, j). In order to apply a likelihood analysis 
we must also specify a noise covariance matrix C^ . We mod- 
elled the noise as a Gaussian random process with variance 
1/n, uncorrelated between pixels, such that Cfj = (1/n) <5y. 

The likelihood of the map, with a particular power spec- 
trum d, is given by 



lnP(C^|x)oc--(x T (C T 



■C N )- 1 x)+Tr[ln{C T + C N )] 



The goal of MLE is to maximize this function, and the 
fastest general method is to use Newton-Raphson iteration 
to find the zeroes of the derivatives in lnP(CV|x) with re- 
spect to C e . We used the MADCAP package (Borrill 1999; 
http://www.nersc.gov/~borrill/cmb/madcap) to derive 
the maximum likelihood banded angular power spectrum 
from the pixelized galaxy map and noise matrix. MADCAP 
is a parallel implementation of the Bond, Jaffe & Knox 
(1998) maximum- likelihood algorithms for the analysis of 
CMB datasets. We ran the analysis software on the super- 
computer Seaborg, administered by the National Energy Re- 
search Scientific Computing Centre (NERSC) at Lawrence 
Berkeley National Laboratory, California. We again applied 
equation 12 to the MADCAP results to correct the measured 
power spectrum for the influence of multiple-component 
sources. 

Boughn & Crittenden (2002) also performed a 
HEALPIX analysis of the NVSS as part of a cross- 
correlation analysis with the CMB searching for evidence of 
the Integrated Sachs- Wolfe effect. From the pixelized map 
they derived an angular correlation function for the NVSS, 
which they used to constrain theoretical models. We com- 
pare their results with ours in Section 6. 



3.5 Testing the methods 

We tested the two methods, direct spherical harmonic and 
maximum likelihood estimation, by generating a dipole dis- 
tribution of N fa 10 s sources over the NVSS geometry, where 
the model dipole possessed the same amplitude and direc- 
tion as that detected in the NVSS (see Blake & Wall 2002b). 
In order to simulate multiple components, we added com- 
panion sources to « 7 per cent of the objects, with the same 
separation distribution as that measured in the NVSS (Blake 
& Wall 2002a). 

The two measurements of the Ce spectrum, corrected 
for multiple components using equation 12, are plotted in 
Figure 4 and are consistent with zero. The Ce measurements 
have been averaged into bands of width A£ = 5, starting 
from 1 = 2. The dipole term I = 1 is of course spuriously 
high, but has a negligible effect on the measured harmon- 
ics at I > 1 - the galaxy dipole (unlike the CMB dipole) 
is only barely detectable in current surveys (Blake & Wall 
2002b). The normalization convention used in Figure 4, and 
the remaining power spectrum plots, is to expand the surface 
overdensity 5 = {a — ao) / ao in terms of spherical harmonics. 
To convert from the definition of Ce in equations 1 and 2 we 
simply divide Ag m by ao (in units of sr _1 ) and hence Ce by 
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Figure 4. The C'e spectrum of a model dipole distribution with 
added multiple-component sources. The spectrum is measured us- 
ing both spherical harmonic analysis (solid circles) and maximum 
likelihood estimation (open circles). The measurements were av- 
eraged into bands of width A£ = 5, starting from 1 = 2. The 
results are consistent with Ce = 0. 



Figure 4 permits a first comparison of the two indepen- 
dent methods of estimating Ce ■ At low £, the variances are in 
excellent agreement. As £ approaches 2 x n s id e , the variance 
of the maximum likelihood method begins to exceed that 
of the spherical harmonic analysis. This occurs as the res- 
olution of the pixelization scheme becomes important: the 
angular pixel size is no longer much less than the character- 
istic angular scale probed by the Ith multipole. The variance 
of the high t bins could be reduced by adopting a finer pix- 
elization, such as n s ido = 64, with the penalty of a rapidly 
increasing requirement of supercomputer time. Given that 
the signal in the NVSS angular power spectrum turns out 
to be confined to I < 40, our optimum pixelization remains 
"-side = 32. Moreover, higher pixel resolution would decrease 
the mean pixel count n in equation 13, rendering a Gaussian 
distribution a poorer approximation for Xi. 



4 RESULTS 

4.1 The angular power spectrum Ce 

Figure 5 plots the NVSS angular power spectrum measured 
for flux-density threshold Si.^ghz = 10 mjy using both 
spherical harmonic analysis and maximum likelihood esti- 
mation. The constant offset due to double sources (equation 
12) has been subtracted and the measurements are averaged 
into bands of width A£ = 5, starting from £mi n = 2. Mea- 
surements are plotted up to £ = 100, although note that the 
variance of the maximum likelihood estimation is increased 
by pixelization effects above £ ~ 50, as discussed in Sec- 
tion 3.5. Table 1 lists the plotted data. Figure 6 displays the 
same data scaled by the usual CMB normalization factor, 
£(£ + 1)/2ty. 

We detect clear signal in the NVSS Ce spectrum at 
multipoles £<A0, well- fitted by a power-law Ce = a£~ b , 
where the best-fitting parameters are a = (2.0 ±0.4) x 10~ 4 , 
b = 1.20 ± 0.06 (fitted to the spherical harmonic analysis 
result) or a = (1.6 ± 0.4) x 10 -4 , b = 1.12 ± 0.08 (fitted to 
the maximum likelihood method result). The amplitude of 




• Spherical harmonics 
° Maximum likelihood 



Multipole I 

Figure 5. The NVSS Ce spectrum for a flux-density threshold 
Si. 4 GHz = 10 m Jyi determined using two different methods. The 
Ce measurements are averaged into bins of width A£ = 5, starting 
from 1 = 2. The solid circles are the result of a spherical harmonic 
analysis and the open circles are derived from a maximum likeli- 
hood analysis. For clarity, the open circles are offset by A£ = — 1. 
The solid line displays the best-fitting power-law to the spherical 
harmonic analysis results. 



Table 1. Data table of banded NVSS Ce values plotted in Figure 
5. The offset due to double sources is ACe = 0.24 x 10~ 5 . 



p ■ 




Spherical harmonic 
C° bs x 10 5 


Maximum likelihood 
C° bs x 10 5 


2 


6 


2.98 ± 1.93 


3.35 ±1.62 


7 


11 


1.46 ± 0.64 


1.19 ±0.64 


12 


16 


0.85 ± 0.44 


0.75 ±0.45 


17 


21 


0.68 ± 0.33 


0.71 ±0.38 


22 


26 


0.79 ± 0.30 


0.41 ±0.31 


27 


31 


0.33 ± 0.23 


0.54 ±0.29 


32 


36 


0.34 ±0.21 


0.49 ± 0.27 


.37 


11 


0.12 ± 0.18 


0.32 ±0.25 


42 


16 


0.18 ± 0.17 


-0.07 ±0.21 


47 


51 


0.30 ±0.17 


0.22 ±0.22 


52 


56 


0.21 ± 0.16 


0.22 ±0.22 


57 


61 


0.26 ± 0.16 


0.31 ±0.23 


62 


66 


0.05 ± 0.14 


0.21 ±0.22 


67 


71 


0.06 ± 0.13 


-0.06 ± 0.21 


72 


76 


0.25 ± 0.14 


0.35 ±0.24 


77 


81 


0.13 ± 0.13 


0.17 ±0.24 


82 


86 


-0.03 ±0.11 


0.01 ±0.24 


87 


91 


-0.05 ±0.11 


-0.15 ± 0.24 


92 


96 


0.08 ± 0.11 


0.09 ± 0.27 


97 


101 


0.12 ± 0.11 


0.11 ±0.27 



the angular power spectrum for radio galaxies is two orders 
of magnitude smaller than that found for optically-selected 
galaxies (e.g. Huterer, Knox & Nichol 2001); this is readily 
explained by the wide redshift range of the NVSS sources, 
which vastly dilutes the clustering signal through the su- 
perposition of unrelated redshift slices. The NVSS signal 
remains ~ 5 orders of magnitude greater than the CMB 
Ce spectrum over the same multipole range, reflecting the 
growth of structure since z = 1100. 

Given the incomplete sky and finite resolution, the mea- 
sured Ce values are not independent. However, it was ar- 
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• Spherical harmonics 
° Maximum likelihood 
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Multipole I 



Figure 6. The same data points and power-law fit as Figure 5, 
scaled by a factor £(£+l) /2n to match the usual CMB convention. 
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Figure 7. The NVSS Ci spectrum measured using spherical har- 
monic analysis for flux-density thresholds of 5 mJy (triangles, 
(7 o = 28.2 deg" 2 ), 10 mJy (solid circles, cr = 16.9 deg" 2 ) and 20 
mJy (open circles, <rn = 9.5 deg -2 ). The 5 mJy points are offset 
by A£ = — 1 and the 20 mJy points are offset by A£ = +1. 



gued in Section 3.2 that the correlations between neigh- 
bouring power spectrum measurements are small. This fact 
was confirmed by the maximum likelihood analysis. The de- 
gree of correlation between neighbouring bins is given by 
the immediately off-diagonal elements of the inverse Fisher 
matrix. This data is generated by the MADCAP software, 
and inspection revealed that the size of the immediately off- 
diagonal matrix elements was ~ 20 times smaller than that 
of the diagonal elements. 

Figure 7 compares the angular power spectra measured 
at flux-density thresholds 5 mJy, 10 mJy and 20 mJy us- 
ing spherical harmonic analysis. The 5 mJy data may be 
affected by systematic surface density gradients. The results 
are consistent with an unchanging underlying power spec- 
trum. This is not surprising; the redshift distribution of ra- 
dio sources does not vary significantly between 5 mJy and 
20 mJy, and the angular correlation function has been found 
not to depend on flux density in this range (Blake & Wall 
2002a). 




|A,j7a 

Figure 8. The distribution of values of \Ai m \ 2 (corrected for par- 
tial sky coverage) for the NVSS for 51 < I < 100. The solid line 
is the prediction for an unclustered galaxy distribution, which 
should be a good approximation given that Ci ~ for this range 
in I. However, the slope is a poor fit due to the influence of 
multiple-component sources. When these are taken into account, 
the observed distribution is accurately reproduced (the dashed 
line). 



4.2 Probability distribution of \Ai m \ 2 

An interesting probe of the galaxy pattern is the distribu- 
tion of values of |A^ m | 2 (see Hauser & Peebles 1973). These 
quantities are measured as part of our spherical harmonic 
analysis (Section 3.2). For a random distribution with sur- 
face density oq over a full sky, the central limit theorem en- 
sures that the real and imaginary parts of Ae m — ^ . Yl m (i) 
are drawn independently from Gaussian distributions such 
that the normalization satisfies \Ai m \ 2 = ctq. It is then easy 
to show that x — \Ai m \ 2 has an exponential probability dis- 
tribution for rn^O: 



P(x) dx ■■ 



exp(— x/a ) 



dx 



(14) 



For a partial sky, l^fml 2 is replaced by \Atm — aoIe m \ 2 / Jt m 
(equation 5). 

Figure 8 plots the distribution of observed values of 
\Ai m — (To I tm\ 2 / Jem- We restrict this plot to the multipole 
range 51 < I < 100: for this range of £, Figure 5 demon- 
strates that Ce ~ and thus the survey is well-described by 
a random distribution with additional multiple components. 
For each i we included the range 1 < m < I (negative values 
of m are not independent). Overplotted on Figure 8 as the 
solid line is the prediction of equation 14. Multiple compo- 
nents cause the slope of the observed exponential distribu- 
tion to be shallower than this prediction. Section 3.3 shows 
that the value of < \Ai m — aoItm\ 2 / Jtm > is increased from 
(Jo to (l + 2e)<7o, where e is the fraction of galaxies split into 
double sources. Thus equation 14 must be amended such 
that P(x) oc exp [— x/(l + 2e)ao]. Assuming that e = 0.07, 
this corrected prediction is plotted on Figure 8 as the dashed 
line and provides a very good fit to the observed distribution. 
This is an independent demonstration that approximately 7 
per cent of NVSS galaxies are split into multiple-component 
sources. Figure 8 also underlines the fact that the imprint 
of clustering on the projected radio sky is very faint. 
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5.1 Theory 



• Spherical harmonics 
° Maximum likelihood 
w(0) mode 




40 60 
Multipole I 

Figure 9. The same data points as Figure 5, overplotted with the 
prediction of equation 3 assuming an angular correlation function 

w(6) = (i x io- 3 )e-°- s . 



4.3 Comparison with w(9) 

The angular correlation function w{9) has been measured 
for the NVSS by Blake & Wall (2002a) and Overzier et al. 
(2003). It is well-described by a power- law w(9) w (1 x 
10 -3 ) 9~ os for angles up to a few degrees. Equation 3 allows 
us to derive the equivalent Ce spectrum if we assume that 
this power-law extends to all angular scales. In Figure 9 we 
overplot the resulting prediction on the measurements and 
find an excellent fit. 

This is initially surprising: the angular correlation func- 
tion is only measurable at small angles 9 < 10°, whereas low 
multipoles £ describe surface fluctuations on rather larger 
angular scales. However, in Section 5 we establish that the 
signal in low multipoles is actually generated at low redshift 
by spatial fluctuations on relatively small scales, similar to 
the spatial scales which produce the signal in w(9). 

Of course the comparison is not entirely straight- 
forward because the two statistics quantify different proper- 
ties of the galaxy distribution. The value of Ce quantifies the 
amplitude of fluctuations on the angular scale correspond- 
ing to £. The value of w(9) is the average of the product of 
the galaxy overdensity at any point with the overdensity at 
a point at angular separation 9 - w(9) depends on angular 
fluctuations on all scales. This is illustrated by the inverse 
of equation 3: 



w{9) 



u £=1 



(2£ + l) C e P e (cost 



For example, a density map constructed using just one mul- 
tipole I possesses a broad angular correlation function. 



5 RELATION TO THE SPATIAL POWER 
SPECTRUM 

The purpose of this Section is to demonstrate that the ob- 
served radio galaxy Ce spectrum can be understood in terms 
of the current fiducial cosmological model. In Section 6 we 
utilize this framework to determine the linear bias factor of 
the radio galaxies, marginalizing over other relevant model 
parameters. 



The angular power spectrum Ce is a projection of the spa- 
tial power spectrum of mass fluctuations at different red- 
shifts, P(k,z), where A; is a co-moving wavenumber. In lin- 
ear perturbation theory, fluctuations with different k evolve 
independently, scaling with redshift according to the growth 
factor D(z). Under this assumption we can simply scale the 
present-day matter power spectrum Po(k) back with red- 
shift: 



P(k,z)=P (k)D(z) 



(15) 



For an O m = 1, Qa = universe, D(z) = (1 + z)^ 1 . For 
general cosmological parameters we can use the approxima- 
tion of Carroll, Press & Turner (1992). If linear theory holds 
then the angular power spectrum Ce can be written in terms 
of the present-day matter power spectrum Po(k) as 



C e = / P (k)W e {k)dk 



(16) 



(see e.g. Huterer, Knox & Nichol 2001, Tegmark et al. 2002) 
where the kernel Wi(k) is given by 



W e (k) = - 



jt(u) f(u/k) du 



(17) 



Here, je is a spherical Bessel function and f(x) is a function 
which depends on the radial distribution of the sources as 



/(*(*)) = 



p(z)D(z)b(z) 
dx/dz 



(18) 



assuming a flat geometry. x(z) is the co-moving radial co- 
ordinate at redshift z, p(z) is the redshift probability distri- 
bution of the sources (normalized such that J p(z) dz = 1) 
and b(z) is a linear bias factor which relates the clustering 
of galaxies to clustering of the underlying mass: 

P ga ,i(k, z) = b(z) 2 P mass (fc, z) 

For the purposes of this analysis we assumed that the linear 
bias does not evolve with epoch and may be represented by 
a constant bias factor b(z) = 60 • 

A useful approximation for the spherical Bessel function 
(which gets better as £ gets larger) is ji(x) « (ir/2£) 1 / 2 5(x — 
£). In this approximation, 



Wt(k) w - f(£/k) 



(19) 



Thus as £ is increased, the kernel just translates along the 
fc-axis. Combining equations 16 and 19 produces the approx- 
imation 



C, 



' ' P (k) f(£/k) 2 dk 



(20) 



or, converting this into an integral over radial co-ordinate, 



Ce 



P (£/x)x~ 2 f(x) 2 dx 



(21) 



In Section 4.1 we found that the measured NVSS Ce spec- 
trum was well-fit by a power-law. Inspecting equations 20 
and 21, a power law for Ce can arise in two ways: 

• If the function f(x) were a power-law f(x) oc x n , then 
equation 20 predicts that Ce oc £ 2 "~ 1 , regardless of the form 
of Po(k). In particular, if f(x) were approximately constant 
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Table 2. Functional form of common clustering statistics for pure 
power-law clustering parameterized by slope 7. 



Statistic 


Dependence 


Validity 


Spatial correlation function 


£(r) oc r ~ 1 


All r 


Angular correlation function 


w(0) oce 1 —! 


Small 6 


Spatial power spectrum 


P(k) OC fcT" 3 


All k 


Angular power spectrum 




High I 



over the relevant scales, then Ce oc £~ , which is a good 
description of the observed NVSS clustering. 

• If the spatial power spectrum were a power-law Po(k) oc 
k", then equation 21 predicts that Ce oc £ n , regardless of the 
form of the radial distribution p(z) (see also Table 2). The 
result Ce oc £ n differs from the case of CMB fluctuations, 
for which n = 1 corresponds to Ce oc l/£(£ + 1) at low 
£ - the well-known Sachs- Wolfe effect. This is because the 
Sachs- Wolfe effect is sensitive to fluctuations in gravitational 
potential, whereas we are probing fluctuations in mass. 

In Section 5.3 we determine that the first of these interpre- 
tations is correct. 



O 




~I0~ J 0.01 0.1 1 

k / h Mpc" 1 

Figure 10. The fiducial model present-day matter power spec- 
trum Po(k) used in this investigation (solid line). This is modified 
from a primordial featureless power-law (dotted line) using the 
transfer function of Eisenstein &; Hu (1998), and from the linear 
power spectrum (dashed line) using the non-linear prescription of 
Peacock & Dodds (1996). The assumed cosmological parameters 
were h = 0.73, U m h 2 = 0.134, n b /Q m = 0.17, n = 0.97 and 
CT 8,lin = !■ 



5.2 Modelling the present-day power spectrum 

We assumed that the primordial matter power spectrum is a 
featureless power-law, P pr i m (A;) = A k n . On very large scales, 
the only alteration to this spectrum in linear theory will be 
an amplitude change due to the growth factor. However, 
during the epoch of radiation domination, growth of fluctu- 
ations on scales less than the horizon scale is suppressed by 
radiation pressure. This process is described by the trans- 
fer function T(k), such that the present-day linear matter 
power spectrum is given by 

Po(k) = P prim (fc) T{k) 2 = Ak n T{kf (22) 

Accurate fitting formulae have been developed for the trans- 
fer function T(k) in terms of the cosmological parameters 
(Eisenstein & Hu 1998), which we employed in our analysis 
(these fitting formulae assume adiabatic perturbations). In 
our fiducial cosmological model, we fixed the values of the 
cosmological parameters at h = i/o/(100 kms -1 Mpc -1 ) = 
0.73, ft m h 2 = 0.134 and Q b /Q, m = 0.17 (Spergel et al. 2003, 
Table 7, column 3). We also chose a primordial spectral in- 
dex n = 0.97 (Spergel et al. 2003), which is close to the 
predictions of standard inflationary models. We consider the 
effect of variations in these parameter values in Section 6. We 
assumed that the Universe is flat, with the remaining energy 
density provided by a cosmological constant S1a = 1 — fi TO . 

There are at least two independent ways of estimating 
the amplitude A in equation 22. Firstly we can use con- 
straints on the number density of massive clusters at low 
redshift, expressed in terms of as, the rms fluctuation of 
mass in spheres of radius R — 8 ft -1 Mpc: 

0% = J Mk) [ji{kR)f dk (23) 

For example, Viana & Liddle (1999) determined the most 
likely value of as using this method to be as = 0.56 Q.^' 47 . 
Alternatively, A can be expressed in terms of the amplitude 



of fluctuations at the Hubble radius, Sh, and constrained by 
measurements of CMB anisotropies on large angular scales: 

A = 2n 2 5 2 H (^) i+n (24) 

For example, Bunn & White (1997) give the best-fitting con- 
straint on 8h and n for flat models (S7 m + Qa = 1) based on 
results of the COBE DMR experiment: 

S H = i.94xlO- 5 fi- - 785 - ' 051nn '" e - a95( ' l - 1) - ' 169( "- 1)2 (25) 

with a maximum la statistical uncertainty of 7 per cent. For 
our fiducial cosmological parameters, a (reasonably) consis- 
tent cosmology is produced if as = 1 (i.e. 5h = 4.7 x 10 ). 
We assumed this fiducial normalization in our model, noting 
that the value of as is in fact degenerate with the amplitude 
of a constant linear bias factor b(z) = bo- 

Equation 22 is only valid ignoring non-linear effects, 
which will boost the value of P(k) on small scales as modes 
commence non-linear collapse. We incorporated non-linear 
corrections using the fitting formula provided by Peacock & 
Dodds (1996). The resulting model power spectrum is shown 
in Figure 10. Strictly, this modification violates the assump- 
tion of linear evolution implicit in equation 15. However, this 
is not significant in our analysis because the small scales for 
which non-linear evolution is important are only significant 
in the projection at z ~ 0. 

5.3 Modelling the radial distribution of NVSS 
sources 

The projection of the spatial power spectrum onto the sky 
depends on the radial distribution of the sources under con- 
sideration, which may be deduced from their redshift distri- 
bution. We only need to know the probability distribution 
p(z) of the sources (equation 18), the absolute normaliza- 
tion is not important. Unfortunately, the radial distribution 
of mjy radio sources is not yet accurately known. The ma- 
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jority of radio galaxies are located at cosmological distances 
(z ~ 1) and their host galaxies are optically very faint. 

However, models exist of the radio luminosity function 
of AGN (i.e. the co-moving space density of objects as a 
function of radio luminosity and redshift), from which the 
redshift distribution at any flux-density threshold can be in- 
ferred. Such models have been published by Dunlop & Pea- 
cock (1990) and Willott et al. (2001). These luminosity func- 
tion models are by necessity constrained by relatively bright 
radio sources (Si.4GHz > 100 mjy) and the extrapolation to 
NVSS flux-density levels must be regarded as very uncer- 
tain. The Willott model is constrained by a larger number 
of spectroscopic redshifts, and the samples of radio sources 
used provide fuller coverage of the luminosity-redshift plane 
(thus the required extrapolation to Si.4GHz = 10 mjy is 
less severe). However, the Willott samples are selected at 
low frequency (151 and 178 MHz), necessitating a large ex- 
trapolation to the NVSS observing frequency of 1.4 GHz. 
The Dunlop & Peacock models are constrained at high fre- 
quencies, but treat steep-spectrum and flat-spectrum radio 
sources as independent populations, which is inconsistent 
with current ideas concerning the unification of radio AGN 
(e.g. Jackson & Wall 1999). In addition, they were computed 
for cosmological parameters Q. m = 1, f?A = rather than 
the currently favoured "ACDM" cosmology. Furthermore, 
none of the aforementioned luminosity function models in- 
corporate starburst galaxies, which contribute in significant 
numbers to the radio galaxy population mix at z<0.1 for 
flux-density threshold 5i.4GHz = 10 mjy. 

Direct measurements of p(z) are currently only achiev- 
able at low redshifts («<0.2), where comparison with 
large optical galaxy redshift surveys is possible (Sadler 
et al. 2002; Magliocchetti et al. 2002). We matched the 
NVSS 10 mjy catalogue with the final data release of the 
2dF Galaxy Redshift Survey (2dFGRS, available online at 
http://msowww.anu.edu.au/2dFGRS/, also see Colless et al. 
2001), in order to estimate p(z) at low redshift. As our fidu- 
cial model, we fitted the resulting redshift histogram with 
the simplest possible function, a constant p(z) — po over the 
range < z < 0.15. In Section 6 we include the effect of vari- 
ations in this model. We only considered the redshift range 
z < 0.15 in this analysis because atz« 0.15, the (very lumi- 
nous) optical counterparts of the 10 mjy NVSS sources be- 
gin slipping below the 2dFGRS magnitude threshold, which 
we verified by plotting the magnitudes of matched 2dFGRS 
galaxies against redshift. Having determined the value of po, 
we created the full redshift distribution by assigning the re- 
maining probability 1 — po over the redshift range z > 0.15 
in proportion to the prediction of the Dunlop & Peacock 
(1990) luminosity function models. For this investigation we 
used the average of the seven models provided by Dunlop & 
Peacock. Assuming the Willott et al. (2001) radial distribu- 
tion for z > 0.15 made a negligible difference to the results 
(because most of the contribution to the Ce spectrum arises 
at low redshifts, see Section 5.4). 

Matching NVSS catalogue entries brighter than 
Si. 4 GHz = 10 mjy with the 2dFGRS database yielded 
Amat = 546 identifications with redshifts z < 0.15, us- 
ing matching tolerance 10 arcsec (see Sadler et al. 2002). 
We restricted the 2dFGRS sample to "high quality" spectra 
(Q > 3). An estimate of the probability of an NVSS source 
being located at z < 0.15 is 



CNVSS 

where ctnvss = 16.9 deg~ 2 is the surface density of NVSS 
sources brighter than 10 mjy, and is the 2dFGRS area 
under consideration, which (owing to the varying angular 
completeness) is not trivial to calculate. We followed Sadler 
et al. (2002) by dividing the number of 2dFGRS galaxies 
contained in the NVSS geometry by the 2dFGRS surface 
density <72dF = 180 deg~ 2 , resulting in an effective area 
A 2dF = 1214 deg 2 . 

The result, po = 0.177 (per unit redshift), is a significant 
underestimate for various reasons: 

• Incompletenesses in the 2dFGRS input catalogue (5%). 

• Input catalogue galaxies unable to be assigned a 2dF 
spectrograph fibre (7%). 

• Observed 2dFGRS spectra with insufficient quality 
(Q < 2) to determine a redshift (8%). 

• Extended radio sources with catalogue entries located 
more than 10 arcsec from the optical counterpart (3%). 

The estimated correction factors in brackets were obtained 
from Colless et al. (2001) and from Carole Jackson (priv. 
comm.). Multiplying these corrections implies a total incom- 
pleteness of 25%, and on this basis we increased the value 
of po to 0.177/0.75 = 0.237. 

Figure 11 plots 2dFGRS-NVSS matches in redshift bins 
of width Az — 0.01, together with the low-redshift fit de- 
scribed above. A constant p(z) is a fairly good approxima- 
tion at low redshifts (z < 0.15). This flat distribution arises 
because the overall redshift distribution is a sum of that 
due to AGN and that due to starburst galaxies; p(z) in- 
creases with z for the AGN, but decreases with z for the 
starbursters. Figure 11 also displays the predictions of the 
luminosity function models of Dunlop & Peacock (1990) and 
Willott (2001). As explained above, the large extrapolations 
involved render these models a poor fit to the redshift distri- 
bution at mjy flux levels, and their use without modification 
at low redshift would have caused significant error. 

5.4 Predicting the Ce spectrum 

We used equation 16 to predict the Ce spectrum from our 
fiducial models of the spatial power spectrum (Section 5.2) 
and the radial distribution of the sources (Section 5.3). We 
found that a good match to the measured angular power 
spectrum resulted if the NVSS sources were assigned a con- 
stant linear bias factor 6o « 1.7 (Figure 12); 6 = 1 provides 
a very poor fit to the results. The bias factor of the radio 
galaxies is analyzed more thoroughly in Section 6. 

We note that these measurements of the radio galaxy Ce 
spectrum at low £ are not directly probing the large-scale, 
small k, region of the power spectrum P(k). Investigation 
of the integrands of equations 20 and 21 revealed that the 
majority of the signal is built up at low redshift, z < 0.1 (see 
Figure 13), where small-scale spatial power is able to con- 
tribute on large angular scales (i.e. contribute to low mul- 
tipoles). Higher redshift objects principally serve to dilute 
the clustering amplitude. This is unfortunate: the potential 
of radio galaxies distributed to z ~ 1 to probe directly the 
large-scale power spectrum is forfeited by projection effects. 
In order to realize this potential, we must measure redshifts 
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Redshift z 

Figure 11. The redshift distribution dN/dz of NVSS sources 
brighter than 10 mjy. The data points are measured values re- 
sulting from cross-matching with the 2dF Galaxy Redshift Survey, 
and the solid line is the best- fitting constant dN/dz = poX c NV ss 
for < z < 0.15. We also plot the predictions of two luminosity- 
function models: Dunlop & Peacock (1990, dashed line) and 
Willott et al. (2001, dotted line). The NVSS-2dFGRS matches 
decline at 2>0.2 due to the falling dN/dz of the optical cata- 
logue. 




• Spherical harmonics 
° Maximum likelihood 
- - Model b =2.0 

Model b =1 .7 

Model b„=1.0 



20 40 
Multipole I 

Figure 12. The predicted Cg spectrum for NVSS sources brighter 
than 10 mjy from a reasonable model spatial power spectrum 
(Section 5.2) and model redshift distribution (Section 5.3). The 
amplitude of the measurements can only be matched if radio 
galaxies are assigned a high bias; model predictions are plotted 
for constant linear bias factors bo = 1 (dotted line), bo = 1-7 
(solid line) and bo = 2 (dashed line). 



for the NVSS sources. A three-dimensional map extending 
to z ~ 1 would directly yield P(k) on large scales, defining 
the "turn-over" sketched in Figure 10. 



6 RADIO GALAXY BIAS FACTOR 

The enhanced radio galaxy bias apparent in Figure 12 is 
consistent with the nature of AGN host galaxies: optically 
luminous ellipticals inhabiting moderate to rich environ- 
ments. Similarly high radio galaxy bias factors have been 
inferred from measurements of the spatial power spectrum 
of low-redshift radio galaxies (Peacock & Dodds 1994) and 
from deprojection of the NVSS angular correlation function 
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Figure 13. Contribution to the model NVSS angular power spec- 
trum as a function of redshift for multipoles £ = 10 (solid line) 
and I = 20 (dashed line), derived from equation 21 using the 
assumptions detailed in Sections 5.2 and 5.3. The Figure plots 
K e (z) where C e = J K e (z) dz. The y— axis in the Figure is scaled 
arbitrarily. The discontinuity at z = 0.15 results from the break 
in the model p(z) at this redshift (see Section 5.3). 



w{0) (Blake & Wall 2002a; Overzier et al. 2003). Further- 
more, Boughn & Crittenden (2002) cross-correlated NVSS 
and CMB overdensities in a search for the Integrated Sachs- 
Wolfe effect (see also Nolta et al. 2003). Their analysis in- 
cluded fitting theoretical models to the NVSS angular cor- 
relation function. They derived good fits for a slightly lower 
linear bias factor than ours, bo = 1.3 — * 1.6. This differ- 
ence is due to the assumed redshift distribution. Boughn & 
Crittenden also used the Dunlop & Peacock (1990) average 
model, but we corrected this model using observational data 
at low redshifts z < 0.15. As can be seen from Figure 11, 
our correction reduces the number of low-redshift sources, 
necessitating a higher bias factor to recover the same angu- 
lar correlations. 

In order to derive a formal confidence interval for the 
linear bias parameter bo we must incorporate the effects 
of uncertainties in the underlying model parameters (by 
marginalizing over those parameters). With this in mind, 
we assumed Gaussian priors for Hubble's constant h — 
0.73 ± 0.03, for the matter density Q. m h 2 = 0.134 ± 0.006, 
and for the primordial spectral index n — 0.97 ± 0.03. The 
widths of these priors were inspired by the cosmological 
parameter analysis combining the WMAP satellite obser- 
vations of the CMB and the 2dFGRS galaxy power spec- 
trum (Spergel et al. 2003, Table 7, column 3) and are a 
good representation of our current knowledge of the cosmo- 
logical model. In addition, we considered variations in the 
model for the radial distribution of NVSS sources at low 
redshift (Section 5.3), using a more general fitting formula 
p(z) = a + bz to describe the probability distribution for 
2 < 0.15. For each pair of values of (a, b) we derived the chi- 
squared statistic between the model and the observations, 
Xp( z ), which we converted into an (unnormalized) probabil- 
ity density P v(z) oc exp(-^ w /2). 

Our model is thus specified by values of (bo, h, Q m h 2 , 
n, a, b) from which we can calculate a model Ct spectrum 
and hence a chi-squared statistic with the observations, Xc e i 
corresponding to a probability density Pc e oc exp (— Xc./S). 
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Figure 14. Normalized probability distribution for bo X <rg, ob- 
tained after marginalizing over the other parameters of the model 
as described in Section 6. The bounds of the 68% confidence re- 
gion are marked by the vertical dashed lines. 



We used the spherical harmonic estimation of the Ci spec- 
trum as the observational data. After multiplying Pc e by 
the redshift distribution probability density Pjv( z ) {a, b) and 
the Gaussian prior probability densities for h, Q m h 2 and n, 
we derived the probability distribution for bo by integrating 
over each of the other parameters. We do not marginalize 
over the normalization of the matter power spectrum, <rg, 
because this quantity is degenerate with bo (using equations 
18, 21 and 23: Ci oc &o<rf). 

The resulting normalized probability distribution for 
6o x 08 is displayed in Figure 14, from which we determined a 
68% confidence region bo as = 1.53 — -> 1.87. When combined 
with the WMAP determination of erg = 0.9 ±0.1 (Spergel 
et al. 2003), we infer that b = 1.89 ± 0.27. 



7 CONCLUSIONS 

This investigation has measured the angular power spec- 
trum of radio galaxies for the first time, yielding consistent 
results through the application of two independent meth- 
ods: direct spherical harmonic analysis and maximum like- 
lihood estimation. The NVSS covers a sufficient fraction of 
sky (~ 80 per cent) that spherical harmonic analysis is very 
effective, with minimal correlations amongst different mul- 
tipoles. The form of the Ci spectrum can be reproduced by 
standard models for the present-day spatial power spectrum 
and for the radial distribution of NVSS sources - provided 
that this latter is modified at low redshift through compari- 
son with optical galaxy redshift surveys. The results strongly 
indicate that radio galaxies possess high bias with respect 
to matter fluctuations. A constant linear bias bo ~ 1.7 per- 
mits a good fit, and by marginalizing over the other pa- 
rameters of the model we deduce a 68% confidence interval 
bo as = 1.53 — > 1.87 where erg describes the normalization of 
the matter power spectrum. We find that the majority of the 
angular power spectrum signal is generated at low redshifts, 
z<0.1. Therefore, in order to exploit the potential of radio 
galaxies to probe spatial fluctuations on the largest scales, 
we require individual redshifts for the NVSS sources. 
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